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Abstract 

Phenomena related to the non-perturbative aspects of strong interactions 
at the LHC are discussed with emphasis on elastic and inelastic soft and hard 
diffraction processes. Predictions for the global characteristics and angular 
distributions in proton-proton collisions with elastic and multiparticle final 
states are given. Potential for discovery of the novel effects related to the 
increasing role of the elastic scattering at the LHC energies and its phys- 
ical implications in diffractive and multiparticle production processes are 
reviewed. 
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Introduction 



The knowledge of hadron structure and hadron interaction dynamics is the ul- 
timate goal of strong interaction theory. Nowadays quantum chromodynamics 
(QCD) is generally accepted as such a theory. Perturbative QCD enables one to 
describe successfully spin-averaged observables at short distances. However, per- 
turbative calculations can not be applied already at the distances larger than 0.1 
fm due to chiral symmetry breaking. Moreover, any real hadron hard interaction 
process involves also long range interaction at some stage. The fundamental prob- 
lems of the strong interactions theory are well known and related to confinement 
and chiral symmetry breaking phenomena. These phenomena have deal with col- 
lective, coherent properties of quarks and gluons. Closely related are the problems 
of the spin structure of a nucleon and helicity non-conservation on the quark and 
hadron levels. 

From the phenomenological point of view all above mentioned phenomena 
take place in the region of diffractive physics. Thus understanding of the diffrac- 
tive interactions plays a fundamental role under the studies of high-energy limit 
of QCD. Studies of diffractive processes are also important in the broad context 
of hadron physics problems |ii6j. It was rather surprising that a high relative prob- 
ability of coherent processes at high energies was revealed in the experiments on 
hard diffraction at CERN 1 1 1 and diffractive events in the deep-inelastic scattering 
at HERA lEllSl. Significant fraction of high-t events among the diffractive events 
in deep-inelastic scattering and in hadron-hadron interactions were also observed 
at HERA [4J and Tevatron [5] respectively. These experimental results renewed 
interest in the further experimental and theoretical studies of diffractive processes 
and stimulated interest to the dedicated QCD experimental studies at the LHC (cf. 
0). 

Among the global problems of strong interactions the total cross-section be- 
haviour and its rise constitute a most important question. There are various ap- 
proaches which provide the total cross-section increase with energy but the rea- 
son leading to such behaviour remains obscure. The nature of the total-cross 
section rising energy dependence currently is not understood since underlying 
microscopic mechanism is dominated by the non-perturbative QCD effects and 
even rather simple question on universality or non-universality of this mechanism 
has no definite answer. An important role here belongs to elastic scattering where 
hadron constituents interact coherently. Single diffraction dissociation is a most 
simple inelastic diffractive process and studies of soft and hard final states in this 
process could become the next step after the elastic scattering. 

Multiparticle production and the global observables such as mean multiplicity 
and its energy dependence alongside with the total, elastic and inelastic cross- 
sections provide a clue to the mechanisms of confinement and hadronization. 
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General principles play an essential role in the nonperturbative sector of QCD 
and unitarity which regulates the relative strength of elastic and inelastic pro- 
cesses is the most important one. Owing to the experimental efforts during recent 
decades it has become evident that coherent elastic scattering process will survive 
at high energies and in particular at the LHC. However, it is not evident: will 
hadron interaction remain to be dominated by multiparticle production? Or at 
some distances elastic scattering channel can become playing a dominating role 
at the LHC energy. This question constitutes an important problem for the back- 
ground estimates for the LHC experiments. We give here arguments in favor of 
the unorthodox point of view, i.e. we would like to discuss possible realization of 
the new scattering mode where elastic scattering prevails at super high energies 
and consider experimental signatures in the studies of hadronic interactions at the 
LHC. Such experiments will be crucial for understanding of microscopic nature 
of the driving mechanism which provide rising cross-section, its possible parton 
structure, high-energy limit of strong interactions and approach to the asymptoti- 
cal region. 

Prevalent role of elastic scattering at very high energies has in some extent 
been implied by the limitations for inelastic processes obtained on the basis of the 
general principles. In particular, it has been shown that the effective interaction 
radius of any inelastic process cannot be greater than the interaction radius of the 
corresponding (i.e. the process with the same particles in the initial state) elastic 
scattering process [8|. 

The appearance of antishadowing would be associated with significant spin 
correlations of the produced particles. We briefly mention some spin related ex- 
perimental possibilities then. 

In general we review here particular problems in hadron interactions which in 
some cases closely or in other cases not that much but related to the phenomena of 
antishadowing. Several original results have already been published in [9J, others 
are discussed here for the first time. 

1 Approach to asymptotical region and 
unitarization methods 

It is always important to know how far the asymptotical region lie. Unfortunately, 
at the moment the answer for the above question can be given in the model de- 
pendent way only and currently there is no an universal criterion. There are many 
model parameterizations for the total cross-sections which use In^ s dependence 
for atot{s)- Such models were used widely since the first CERN ISR results had 
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appeared'. This implies the saturation of the Froissart-Martin bound, however 
with coefficient in front of In^ s which is lower than the asymptotical bound ^. On 
the other side the power-like parameterizations of atotis) disregard the Froissart- 
Martin bound and considers it as a matter of the very distant asymptopia. Both 
approaches provide successful fits to the experimental data in the available energy 
range and even lead to similar predictions for the LHC energies. 

However, it is not clear whether the power-like energy dependence would 
obey unitarity bound for the partial-wave amplitudes at the LHC energies and be- 
yond. Meanwhile as it is mentioned unitarity is an important principle which is 
needed to be fulfilled anyway. The most straightforward way is to construct an 
amplitude which ab initio satisfy unitarity. But the most common way consists 
in use of an unitarization procedure of some input power-like "amplitude". Uni- 
tarization provides a complicated energy dependence of atot{s) which can be ap- 
proximated by the various functional forms depending on particular energy range 
under consideration. These forms providing a good description of the experimen- 
tal data in the limited energy range have nothing to do with the true asymptotical 
dependence In^ s. Of course, unitarization will lead to the In^ s dependence but 
only at s — s> oo. 

Unitarity of the scattering matrix SS~^ = 1 implies, in principle, an existence 
at high energies s > sq, where sq is some threshold, of the new scattering mode 
- antishadow one. It has been revealed in [11 J and described in some detail (cf. 
lfT2ll and references therein) and the most important feature of this mode is the 
self-damping of the contribution from the inelastic channels. 

We argue here that the antishadow scattering mode could be definitely revealed 
at the LHC energies and give quantitative and qualitative predictions based on the 
rational unitarization, i.e. U -matrix unitarization method [ 13 1. There is no univer- 
sal, generally accepted method to implement unitarity in high energy scattering 
and as a result of this fact a related problem of the absorptive corrections role 
and their sign has a long history (cf. lfT4ll and references therein). However, a 
choice of particular unitarization scheme is not just a matter of taste. Long time 
ago the arguments based on analytical properties of the scattering amplitude were 
put forward ifTSi in favor of the rational form of unitarization. It was shown that 
this form of unitarization reproduced correct analytical properties of the scattering 
amplitude in the complex energy plane much easier compared to the exponential 
form, (simple eikonal singularities would lead to an essential singularities in the 
amplitude). In potential scattering the eikonal (exponential) and tZ-matrix (ra- 
tional) forms of unitarization correspond to two different approximations of the 

'The first model which provide In^ s dependence for total cross-section was developed by 
Heisenberg 1101 . 

^It is known but not often mentioned fact, that the amplitude which provides an exact saturation 
of the Froissart-Martin bound does correspond to pure elastic scattering. 
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scattering wave function, which satisfy the Schrodinger equation to the same or- 
der. Rational form of unitarization corresponds to an approximate wave function 
which changes both the phase and amplitude of the wave. This form follows from 
dispersion theory. It can be rewritten in the exponential form but with completely 
different resultant phase function, and relation of the two phase functions is given 
in [il5J. 



2 Unitarity: particle production and 
elastic scattering 

In the impact parameter representation the unitarity relation written for the elastic 
scattering amplitude /(s, b) at high energies has the form 

Imf{s,b) = \f{s,b)\^ + r]{s,b) (1) 

where the inelastic overlap function r](s,b) is the sum of all inelastic channel con- 
tributions. It can be expressed as a sum of n-particle production cross-sections at 
the given impact parameter 

Vis,b) = ^anis,b). (2) 

n 

The impact parameter b has a simple geometrical meaning as the distance in the 
transverse plane between the centers of the two colliding hadrons. Unitarity equa- 
tion has the two solutions for the case of pure imaginary amplitude: 

f(s,b)='-[l-y/l-Ms,b)], (3) 



f{s,b)='-[l + y/l-4r^{sM (4) 

Almost everywhere the second solution is not taken into account, since f{s, b) 
and r]{s,b) — > Oatfe — > oo. However, there is nothing wrong with the sec- 
ond solution in the limited region of impact parameters. Existence of the second 
solution leads to interesting experimental predictions and should be taken into 
account. Both solutions of unitarity are naturally reproduced by the rational (f/- 
matrix) form of unitarization. In the [/-matrix approach the form of the elastic 
scattering amplitude in the impact parameter representation is the following: 

1 — iU{s, b) 
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U{s, h) is the generalized reaction matrix, which is considered to be an input dy- 
namical quantity similar to eikonal function. It is worth noting that transition to 
antishadowing at small impact parameters can be incorporated into eikonal uni- 
tarization, however, the latter should have a very peculiar form. Inelastic overlap 
function is connected with t/(s, 6) by the relation 

Construction of the particular models in the framework of the U -matrix ap- 
proach proceeds the common steps, i.e. the basic dynamics as well as the notions 
on hadron structure being used to obtain a particular form for the [/-matrix. U- 
matrix unitarization scheme and eikonal scheme lead to different predictions for 
the asymptotical behaviour of inelastic cross-section and for the ratio of elastic 
to total cross-section. This ratio in the [/-matrix unitarization scheme reaches its 
maximal possible value at s — > oo, i.e. growth of the elastic cross-section is to 
be steeper than the growth of the inelastic cross-section beyond some threshold 
energy 

which reflects in fact that the upper bound for the partial-wave amplitude in the 
[/-matrix approach (unitarity limit) is |/;| < 1 while the bound for the case of 
imaginary eikonal is (black disk limit): |/;| < 1/2. When the amplitude exceeds 
the black disk limit (in central collisions at high energies) then the scattering at 
such impact parameters turns out to be of an antishadow nature. In this antishadow 
scattering mode the elastic amplitude increases with decrease of the inelastic chan- 
nels contribution. 

It is worth noting that the shadow scattering mode is considered usually as 
the only possible one. But as it was already mentioned existence of the second 
solution of unitarity in the limited range of impact parameters is completely lawful 
and an antishadow scattering mode should not be excluded. Antishadowing can 
occur in the limited region of impact parameters h < R{s) (while at large impact 
parameters only shadow scattering mode can be realized. Shadow scattering mode 
can exist without antishadowing, but the opposite statement is not valid. 

Appearance of the antishadow scattering mode is consistent with the basic idea 
that the particle production is the driving force for elastic scattering. Indeed, the 
imaginary part of the generalized reaction matrix is the sum of inelastic channel 
contributions: 

lmU{s,b)=J2Un{s,b), (8) 

n 

where n runs over all inelastic states and 

Un{s,b) = J dK\Un{s,b,{CnW (9) 
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and dVn is the n-particle element of the phase space volume. The functions 
Un{s,b, {Cn}) are determined by the dynamics of 2 ^ n processes. Thus, the 
quantity ImU (s, b) itself is a shadow of the inelastic processes. However, unitar- 
ity leads to self-damping of the inelastic channels [TS] and increase of the func- 
tion Imf/(s, b) results in decrease of the inelastic overlap function 77(5, b) when 
ImU {s,b) exceeds unity. 

Respective inclusive cross-section ifTTlfTSl which takes into account unitarity 
in the direct channel has the form 

— = Svr / bdb- — — -,, 



vr / bdb- — \;,; [,,^ . (10) 



The function I(s, b, C) inEq. dlOb is expressed via the functions f/„(s, b, C, {Cn-i}) 
determined by the dynamics of the processes hi + h2 ^ + Xn-i'. 

I{s,b,C) = J2n I rfr„|f/„(s,6,C,{Cn-i})|' (11) 



and 

j I{s, b, C)dC = n{s, b)lmU{s, b). (12) 

The kinematical variables C, (x and pj_, for example) refer to the produced 
particle h-^ and the set of variables {Cn-i} describe the system of n — 1 
particles. 

Let us consider now transition to the antishadow scattering mode, which was 
revealed in flllll . With conventional parameterizations of the [/-matrix (which 
provide rising cross-sections) inelastic overlap function increases with energies 
at modest values of s. It reaches its maximum value r]{s, = 0) = 1/4 at some 
energy s = Sq and beyond this energy the antishadow scattering starts to de- 
velop at small values of b first. The region of energies and impact parameters 
corresponding to the antishadow scattering mode is determined by the conditions 
Im/(s, b) > 1/2 and 77(5, b) < 1/4. The quantitative analysis of the experimen- 
tal data |fT9ll gives the threshold value of energy: ^/sq ~ 2 TeV. This value is 
confirmed by the another model considerations ll20ll . 

Thus, the function r]{s, b) becomes peripheral when energy is increasing be- 
yond s = Sq. At such energies the inelastic overlap function reaches its maximum 
value at 6 = R{s) where R{s) is the interaction radius. So, beyond the transition 
energy there are two regions in impact parameter space: the central region of an- 
tishadow scattering at 6 < R{s) and the peripheral region of shadow scattering at 
b > R{s). The impact parameter dependence of the inelastic channel contribution 
77(5, 6) at s > So are represented in Fig. 1 for the case of standard unitarization 
scheme and for the unitarization scheme with ani shadowing. 
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Figure 1 : Impact parameter dependence of the inelastic overlap function in the 
standard unitarization scheme (left panel) and in the unitarization scheme with 
antishadowing (right panel). 



The region of the LHC energies is the one where antishadow scattering mode is 
to be presented. It will be demonstrated in the next section that this mode can be 
revealed directly measuring aei{s) and <Jtot{s) and not only through the analysis 
in impact parameter representation. 



3 Approach to asymptotics in the ^/-matrix model 

To get the numerical estimates we shall use the following ansatz for the general- 
ized reaction matrix 

N 

exp(-M6/0 = igi-s) exp(-M6/0, (13) 

where M = "^^^i rag. Here mg is the mass of constituent quark, which is taken 
to be 0.35 GeV, N is the total number of valence quarks in the colliding hadrons, 
i.e. iV = 6 for pp-scattering. The value for the other parameters were obtained 
in 113 and have the following values g = 0.24, ^ = 2.5, a = 0.56 ■ 10"^. 
With these values of parameters the model provides satisfactory description of the 
available experimental data for the forward elastic pp-scattering. To obtain the 
above explicit form for the function U{s, h) we used chiral quark model for U- 
matrix [21J where U (s, h) is chosen as a product of the averaged quark amplitudes 

N 

U{s,h) = \[{fQ{sM) (14) 

Q=l 



U{s,h) = ig 



1 + a 
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in accordance with the assumed quasi-independence of valence quark scattering. 
The 6-dependence of the function (/g) has a simple form (/g) oc exp(— mgfe/i^) 
which correspond to quark interaction radius rq = ^/mg. 
For the LHC energy ^/s = 14 TeV we have 

atot ^ 230 mb (15) 

and 

a,i/atot ^ 0.67. (16) 
Thus, the antishadow scattering mode could be discovered at LHC by measuring 




Figure 2: Total and ratio of elastic to total cross-sections of pp-interactions 



(J el /cF tot ratio which is greater than the black disc value 1/2 (cf. Fig. |2l). 

However, the LHC energy is not in the asymptotic region yet; the total, elastic 
and inelastic cross-sections behave like 



atot,ei oc In 



Asymptotical behavior 



g\l + a 



N 



atnei OC In 



g\l + a 



mq 



N 



atot,ei OC In s, ainei oc In s 



(17) 



(18) 



is expected at > 100 TeV. 

Another predictions of the chiral quark model is decreasing energy depen- 
dence of the the cross-section of the inelastic diffraction at s > sq. Decrease of 
diffractive production cross-section at high energies (s > sq) is due to the fact 
that r]{s, b) becomes peripheral at s > sq and the whole picture corresponds to the 
antishadow scattering at 6 < R(s) and to the shadow scattering atb > R{s) where 
R(s) is the interaction radius: 



da, 



diff 



Ml 



(19) 
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The parameter g* < 1 is the probability of the excitation of a constituent quark 
during interaction. Diffractive production cross-section has the familiar 
dependence which is related in this model to the geometrical size of excited con- 
stituent quark. 

At the LHC energy ^/s = 14 GeV the single diffractive inelastic cross- 
sections is limited by the value adiff{s) < 2.4 mb. 

The above predicted values for the global characteristics of pp - interactions at 
LHC differ from the most common predictions of the other models. First, the total 
cross-section is predicted to be twice as much of the common predictions in the 
range 95-120 mb [22J and it even overshoots the existing cosmic ray data. How- 
ever, extracting total proton-proton cross sections from cosmic ray experiments 
is model dependent and far from straightforward (see, e.g. ll23l and references 
therein). Those experiments measure the attenuation lengths of the showers ini- 
tiated by the cosmic particles in the atmosphere and are sensitive to the model 
dependent parameter called inelasticity. They do not provide any information on 
elastic scattering channel. 

4 Angular structure of elastic scattering and diffrac- 
tion dissociation 

Elastic scattering amplitude F(s,t) is determined by the singularities in the im- 
pact parameter complex plane P = b'^ plane. It has poles which positions are 
determined by the solutions of the following equation: 

1 + U{s,p)=0 (20) 

and the branching point at /3 = (cf. Fig. 3), i.e. 

F{s,t)=Fp{s,t)+F,{s,t). 



Contribution of the poles located at the points 



R(s) + i^vrn 



n 



±1,±3,... 



where 



R{s) 



M 



In g{s) 



determines the elastic amplitude in the region \t\/s <^ 1 (t ^ 0). The amplitude 
in this region can be represented in a form of series over the parameter r(-\/^): 



(21) 



k=l 
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Figure 3: Singularities of scattering amplitude in the complex /9-plane 
where the parameter r decreases exponentially with -y/^: 



This series provides diffraction peak and dip-bump structure of the differential 
cross-section in elastic scattering. In the region of moderate t it is sufficient to 
keep few or even one of the terms of series Eq. (|2T1) . The differential cross-section 
in this region has well known Orear behavior. For elastic scattering amplitude 
F{s, t) the pole and cut contributions are decoupled dynamically when g(s) oo 
at s — s> oo. At large momentum transfer (s — oo, \t\/s - fixed) the contribution 
from the branching point (P = 0) is a dominating one. The angular distribution in 
this region has the power dependence 



There is a direct interrelation of the power law behavior of the differential cross 
sections of large angle scattering with rising behavior of the total cross sections at 
high energies. 

Similarity between elastic and inelastic diffraction in the t-channel approach 
suggests that the latter one would have similar to elastic scattering behavior of the 
differential cross-section. However, it cannot be taken for granted and e.g. trans- 
verse momentum distribution of diffractive events in the deep-inelastic scattering 
at HERA shows a power-like behavior without apparent dips ll24l . Similar behav- 
ior was observed also in the hadronic diffraction dissociation process at CERN 
Q where also no dip and bump structure was observed. Angular dependence of 
diffraction dissociation together with the measurements of the differential cross- 
section in elastic scattering would allow to determine the geometrical properties 





(22) 
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of elastic and inelastic diffraction, their similar and distinctive features and origin. 
In the U -matrix approach the impact parameter amplitude of diffraction dissocia- 
tion Fdiff{s, b, Mx) can be written in the pure imaginary case as a square root of 
the cross-section, i.e. 

Fdiffis, b, Mx) = sJUdiff{s,b,Mx)/[l + U{s, b)] (23) 
and the amplitude Fdiff{s,t, Mx) is 

Fd,ff{s,t,Mx) = ^ bdbJo{bV^)^Ud,ff{s, b, Mx)/[l + U{s, b)]. (24) 

The corresponding amplitude Fdiff{s,t,Mx) can be calculated analytically. To 
do so we continue the amplitudes Fdiff{s,P, Mx), P = to the complex 
plane and then Fdiff{s,t, Mx) can be represented as a sum of the pole contribu- 
tion and the contribution of the cut |25|: 

Fdiff{s,t,Mx) = Fdiff,pis,t,Mx) + FdiffAs,t,Mx) (25) 

The situation is different in the case of diffraction production. Instead of dy- 
namical separation of the pole and cut contribution discussed above we have a 
suppression of the pole contribution at high energies since at fixed t 

F'd.ff,p = 0{s[g{s)]~^ ln'/'g{s)), F,,;;,e = 0{s[g{s)]-^). (26) 
Therefore, at all t values we will have 

Fdiff{s,t,Mx) ^ Fdiff,cis,t,Mx), (27) 

where 

Fd,ff,c{s,t,Mx) ^ tg*g-'/\s){l - -^)-^'\ (28) 

where Mx = {Mx — M — l)/2^. This means that the differential cross-section 
of the diffraction production will have smooth dependence on t with no apparent 
dips and bumps 

^""^'ff oc fl - ^)-3 (29) 
dtdMl^^^ Ml' • 

It is interesting to note that at large values of Mx ^ M the normalized differen- 



1 do- 
scaling behavior 



tial cross-section ^ dUM^ is the value of cross-section at t = 0) will exhibit 



11 



Figure 4: Scaling behavior of the normalized differential cross-section 



1 da 



(TO dtdhP ■ 



and explicit form of the function /( 



t/Mj.) is the following 



f{-t/Ml) 



(1 - Aih/Ml)-^. 



(31) 



This dependence is depicted on Fig. 4. 

The above scaling has been obtained in the model approach, however it might 
have a more general meaning. 

The angular structure of diffraction dissociation processes given by Eq. (l29t 
takes place at high energies where g{s) > 1 while at moderate and low energies 
where g{s) < 1 the both contributions from poles and cut are significant. In this 
region 



where the parameter r(v^^) is the same as it is in the elastic scattering. Thus at 
low energies the situation is similar to the elastic scattering, i.e. diffraction cone 
and possible dip-bump structure should be present in the region of small values of 
t, and the overall behaviour of the differential cross-section will be rather compli- 
cated and incorporates diffraction cone, Orear type and power-like dependencies. 

However, at high energies a simple power-like dependence on t Eq. (l29b is 
predicted. It was shown that the normalized differential cross-section has a scaling 
form and only depends on the ratio —t/Mx at large values of M|. 

In fact, our particular comparative analysis of the poles and cut contributions 
has very little with the model form of the fZ-matrix. This is why it has a more 
general meaning. 

At the LHC energies the diffractive events with the masses as large as 3 TeV 
could be studied. It would be interesting to check this prediction at the LHC where 



oo 



(32) 



k=l 
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the scaling and simple power-like behavior of diffraction dissociation differential 
cross-section should be observed. Observation of such behavior would confirm 
the diffraction mechanism based on excitation of the complex hadronlike object - 
the constituent quark. 

5 Angular distributions of leading protons in cen- 
tral production processes 

It was proposed in f26l to study influence of centrally produced particles on the 
angular distribution of leading protons in the processes with two rapidity gaps, 
which are also known as a double Pomeron exchange (dpe) processes: 

p + p^p + X+p. (33) 

In what follows we will consider symmetrical case h = t2 = t. It is interesting to 
study how the diffractive pattern observed in elastic scattering will be changed in 
the corresponding processes with centrally produced particles. 

First of all elastic scattering and process (l33b are very different from the point 
of view of unitarity equation The process (l33t is one of the many contributing 
processes to the function r](s, h) and consequently into the amplitude of elastic 
scattering /(s, h). In the [/-matrix approach the impact parameter amplitude of 
the process (l33t Fdpds, b, Q can be written in the pure imaginary case according 
to JTOt as a square root of the cross-section, i.e. 

Fdpeis, b, C) = VHs,b,C)/[l + U{s, b)] (34) 
and the amplitude Fdpe{s, t,C) is 

?s r°° , , 

F,p,{s,t,0 = - bdbUb^ri)^I{s, 6, C)/[l + U{s, b)]. (35) 
^ Jo 

The variable t is the momentum transfer to one of the protons, while the variable 
C, is related to the system of particles X or to one particle from this system. Using 
relation (fT2b we can represent /(s, 6, C) in the form 

/(s,6,C) = $(s,6,C)Im[/(s,6), (36) 

where 

j ^{s,bX)dC = n{s,b). (37) 

For the mean multiplicity we suppose that the multiplicity of the centrally pro- 
duced particles is given by the following expression 

n{s,b)=(3N^{s)Dc{h), (38) 
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where the function Dc{b) describes distribution of two hadron condensates clouds 
in the overlapping region. Arguments in favor of such form and assumed hadron 
structure are desribed in the next section. The corresponding amplitude Fdpe{s, t, () 
can be calculated analytically and calculations are similar to the case of elastic 
scattering amplitude and amplitude of diffraction dissociation. It is necessary to 
continue the amplitudes Fd-peis, P, () {(3 = b^), to the complex /5-plane and trans- 
form the Fourier-Bessel integral over 

impact parameter into the integral in the complex P - plane over the contour C 
which goes around the positive semiaxis. The amplitude Fdpe{s, P,C) has the 
poles and a branching point at /5 = 0. Therefore it can be represented as a sum of 
the pole contributions and the contribution of the cut: 

Pdpe Fdpe,p (s,t,C)+Fdpe,c(s,t,C) (39) 

Then, using relation (l37t and assuming 

$(5,6,0 =n(s, 6)0(0 

we obtain that 

Fdpe,p--s Yl Hs,^/K)(p{0]^/KKo{^/tK), (40) 

n=±l,±3,... 

i.e. diffractive pattern of leading protons will depend on the distribution of mean 
multiplicity in impact parameter of the centrally produced particles. Using for 
the mean multiplicity results described in Section 6, the amplitude (l40b can be 
rewritten in the form 

Fdpe,p ss~'^''^U{0 Yl e^™^V^i^o(V^) (41) 

n=±l,±3,... 

Thus, the presence of oscillating factor e*'^" would lead to significant differ- 
ences in the diffractive patterns of leading protons in the processes and elastic 
scattering. Indeed, at small values of t all terms of the series (l4TT) are important. 
In elastic scattering the summation over all n leads to the exponential behavior of 
the differential cross section ||2T1 : 

^ oc exp (B(s)t) , B(s) oc In^ s. (42) 
at 

However, for the process (l33t the terms with the large values of n will be sup- 
pressed due the oscillation factor. Thus, we expect that the Orear type behaviour 
will lake place already at low values of t and differential cross section would have 
the following t-dependence already at small and moderate values of t: 

oc exp — T^\r^ ■ (43) 



dtdi \ M 
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6 Multiparticle production and antishadowing 



The region of the LHC energies is the one where new, antishadow scattering mode 
can be observed. Immediate question arises on consistency of the antishadowing 
with the growth of mean multiplicity in hadronic collisions with energy. More- 
over, many models and the experimental data suggest a power-like energy depen- 
dence of mean multiplicity'' and a priori the compatibility of such dependence 
with antishadowing is not evident. 

Now we turn to the mean multiplicity and consider first the corresponding 
quantity in the impact parameter representation. As it follows from (fTTT) and (fT2t 
the n-particle production cross-section cr„(s, b) 



Then the probability 



P,.(..6).^^ = ^%^. (45) 

ainei{s,b) lmU{s,b) 



Thus, we observe the cancellation of unitarity corrections in the ratio of the 
cross-sections (t„(s, b) and aineiis, b). Therefore the mean multiplicity in the im- 
pact parameter representation 



n 



is not affected by unitarity corrections and therefore cannot be proportional to 
r]{s, b). This conclusion is consistent with Eq. (fT2b . The above mentioned propor- 
tionality is a rather natural assumption in the framework of the geometrical mod- 
els, but it is in conflict with the unitarity. Because of that the results [30 1 based 
on such assumption should be taken with precaution. However, the above cancel- 
lation of unitarity corrections does not take place for the quantity n{s) which we 
address now. 

We use a model for the hadron scattering described in ||2ll which is based 
on the ideas of chiral quark models. The picture of a hadron consisting of con- 
stituent quarks embedded into quark condensate implies that overlapping and in- 
teraction of peripheral clouds occur at the first stage of hadron interaction (Fig. 
121). Nonlinear field couplings could transform then the kinetic energy to the inter- 
nal energy and mechanism of such transformation was discussed by Heisenberg 



^Recent discussions of the problems of multiparticle production processes and rising mean 
hadronic multiplicity dependence can be found in 1.27.1 and L28J 
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ifTOll and Carruthers OTl . As a result massive virtual quarks appear in the over- 
lapping region and some effective field is generated. Valence constituent quarks 
located in the central part of hadrons are supposed to scatter simultaneously in a 
quasi-independent way by this effective field. 




Figure 5: Schematic view of initial stage of the hadron interaction. 

Massive virtual quarks play a role of scatterers for the valence quarks in elastic 
scattering and their hadronization leads to production of secondary particles in the 
central region. To estimate number of such scatterers one could assume that part 
of hadron energy carried by the outer condensate clouds is being released in the 
overlap region to generate massive quarks. Then this number can be estimated by: 

iV(s, b) oc ~ (^q))v^ jjh, ^ j^h, ^ N^^s)Dc{b), (46) 

where mg - constituent quark mass, (fcg) - average fraction of hadron energy car- 
ried by the constituent valence quarks. Function D'^ describes condensate distri- 
bution inside the hadron h, and b is an impact parameter of the colliding hadrons. 

Thus, N{s, b) quarks appear in addition to = + rih^ valence quarks. 
In elastic scattering those quarks are transient ones: they are transformed back 
into the condensates of the final hadrons. Calculation of elastic scattering am- 
plitude has been performed in lIlTl . However, valence quarks can excite a part 
of the cloud of the virtual massive quarks and these virtual massive quarks will 
subsequently fragment into the multiparticle final states. Such mechanism is re- 
sponsible for the particle production in the fragmentation region and should lead 
to strong correlations between secondary particles. It means that correlations exist 
between particles from the same (short-range correlations) and different clusters 
(long-range correlations) and, in particular, the forward-backward multiplicity 
correlations should be observed. This mechanism can be called as a correlated 
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cluster production mechanism. Evidently, similar mechanism should be signif- 
icantly reduced in e"^e^-annihilation processes and therefore large correlations 
are not to be expected there. 

As it was already mentioned simple (not induced by interactions with valence 
quarks) hadronization of massive N{s,b) quarks leads to formation of the multi- 
particle final states, i.e. production of the secondary particles in the central region. 
The latter should not provide any correlations in the multiplicity distribution. 

Remarkably, existence of the massive quark-antiquark matter in the stage pre- 
ceding hadronization seems to be supported by the experimental data obtained at 
CERN SPS and RHIC (see [321 and references therein). 

Since the quarks are constituent, it is natural to expect direct proportionality 
between a secondary particles multiplicity and number of virtual massive quarks 
appeared (due to both mechanisms of multiparticle production) in collision of the 
hadrons with given impact parameter: 

n{s, h) = a{nn, + n;,JiVo(s)D^(6) + l3No{s)Dc{h), (47) 

with constant factors a and (3 and 

Dpih) = Dq® Dc, 

where the function Dq (b) is the probability amplitude of the interaction of valence 
quark with the excitation of the effective field, which is in fact related to the quark 
matter distribution in this hadron-like object called the valence constituent quark 
ll21il . The mean multiplicity n(s) can be calculated according to the formula 

n{s) = /o°° ^(^' ^)^(^' ^)^^^ (48) 
Jo° vi^^ h)hdh 

It is evident from Eq. (l48b and Fig. 1 that the antishadow mode with the peripheral 
profile of 77(5, h) suppresses the region of small impact parameters, and the main 
contribution to the mean multiplicity is due to peripheral region of 6 ~ R{s). 

To make explicit calculations we model for simplicity the condensate distri- 
bution Dc{h) and the impact parameter dependence of the probability amplitude 
Dqip) of the interaction of valence quark with the excitation of the effective field 
by the exponential forms, and thus we use exponential dependencies for the func- 
tions Dp{h) and Dc{h) with the different radii. Then the mean multiplicity 

h) = aNois) expi-b/Rp) + pNo{s) exp{-b/Rc). (49) 



After calculation of the integrals (I48t we arrive to the power-like dependence 
of the mean multiplicity n(s) at high energies 

n{s) = as^^ + bs^'^ , (50) 
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where 

There are four free parameters in the model, «, (3 and Rp, Re, and the freedom 




101 10^ 103 104 



s'" (GeV) 

Figure 6: Energy dependence of mean multiplicity, theoretical curve is given by 
the equation n{s) = as^ (a = 2.328, S = 0.201); experimental data from the Refs. 

01. 

in their choice is translated to a, b and Sp, Sc- The value of ^ = 2 is fixed from 
the data on angular distributions lIlTIl and for the mass of constituent quark the 
standard value mg = 0.35 GeV was taken. However, fit to experimental data 
on the mean multiplicity leads to approximate equality Sp — Sc and actually 
Eq. (l50b is reduced to the two-parametric power-like energy dependence of mean 
multiplicity 

n = as^ , 

which is in good agreement with the experimental data (Fig. |6l). Equality Sp — Sc 
means that variation of the correlation strength with energy is weaker than the 
power-like one and could be described, e.g. by a logarithmic function of energy. 
From the comparison with the data on mean multiplicity we obtain that 5 ~ 0.2, 
which corresponds to the effective masses, which are determined by the respective 
radii (M = 1/R), Mc ^ Mp ^ 0.3mQ, i.e. Mp Mc ^ m^. 

The value of mean multiplicity expected at the LHC energy (y/s = 14 TeV) 
is about 110. It is not surprising that it is impossible to differentiate contributions 
from the two mechanisms of particle production at the level of mean multiplicity. 
The studies of correlations are necessary for that purpose. 

Multiplicity distribution P„(s, 6) and mean multiplicity n(s, 6) in the impact 
parameter representation have no absorptive corrections, but since antishadow- 
ing leads to suppression of particle production at small impact parameters and 
the main contribution to the integral multiplicity n(s) comes from the region of 
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b ~ R{s). Of course, this prediction is to be valid for the energy range where 
antishadow scattering mode starts to develop and is therefore consistent with the 
"centrality" dependence of the mean multiplicity observed at RHIC ll35l . 

It would be interesting to note that due to peripheral form of the inelastic over- 
lap function the secondary particles will be mainly produced at impact parameters 
b ~ R{s) and thus will carry out large orbital angular momentum 

L~i?(s)^n(s,6 = /?(s)). 

To compensate this orbital momentum spins of secondary particles should become 
lined up, i.e. the spins of the produced particles should demonstrate significant 
correlations when the antishadow scattering mode appears. Similar conclusion on 
spin correlations of secondary particles has been made long time ago by Yang and 
Chou under utilization of the concept of partition temperature in the geometric 
model for hadron production ll36l . 

It is also worth noting that no limitations follow from the general principles for 
the mean multiplicity, besides the well known one based on the energy conserva- 
tion law. Having in mind relation (l49t . we could say that the obtained power-like 
dependence which takes into account unitarity effects could be considered as a 
kind of a saturated upper bound for the mean multiplicity growth. 

Elastic scattering domination at the LHC and the appearance of the antishadow 
scattering mode implies a somewhat unusual scattering picture. At high energies 
the proton should be represented as a very loosely bounded composite system 
and it appears that this system has a high probability to reinstate itself only in 
the central collisions where all of its parts participate in the coherent interactions. 
Therefore the central collisions are responsible for elastic processes while the pe- 
ripheral ones where only few parts of weekly bounded protons are involved result 
in the production of the secondary particles. This leads to the peripheral impact 
parameter profile of the inelastic overlap function. Such evolution could be ac- 
complished with spin correlations of the produced particles. 

7 Polarization measurements 

In soft hadronic interactions significant single- spin effects could be expected since 
the helicity conservation does not work for interactions at large distances, once 
the chiral SU (3)^ x SU{3) r symmetry of the QCD Lagrangian is spontaneously 
broken in the real world. Thus, studies of the p^-dependence of the one-spin 
asymmetries can be used as a way to reveal a transition from a non-perturbative 
phase (P 7^ 0) to the perturbative one (P = 0). The essential point here is an 
assumption that at short distances the vacuum is perturbative. However, the very 



19 



existence of the above transition can not be taken for granted since the vacuum, 
even at short distances, could be filled up with the fluctuations of gluon and quark 
fields. The measurements of the one-spin transverse asymmetries and polarization 
is an important probe of the chiral structure of the effective QCD Lagrangian. 

At the same time we can note that polarization effects as well as some other 
recent experimental data demonstrate that hadron interactions have a significant 
degree of coherence. Experimentally, spin asymmetries increase at high transverse 
momentum in elastic scattering and are flat in inclusive processes. 

It is interesting to note that on the base of the model in ll37ll one should ex- 
pect a zero polarization in the region where quark-gluon plasma (QGP) has been 
formed, since chiral symmetry is restored and there is no room for quasiparti- 
cles such as constituent quarks. Thus, the absence or strong diminishing e. g. 
of transverse hyperon polarization can be used as a signal of QGP formation in 
heavy-ion collisions. This prediction should also be valid for the models based 
on confinement, e.g. the Lund and Thomas precession model. We could use a 
vanishing polarization of e. g. A-hyperons in heavy ion collisions as a sole result 
of QGP formation provided the corresponding observable is non-zero in proton- 
proton collisions. The prediction based on this observation would be a decreasing 
behavior of polarization of A with the impact parameter in heavy-ion collisions in 
the region of energies and densities where QGP was produced: 

Pa(6)^0 at b^O, (51) 

since the overlap is maximal at 6 = 0. The value of the impact parameter can be 
controlled by the centrality in heavy-ion collisions. The experimental program 
could therefore include measurements of A-polarization in pp-interactions first, 
and then if a significant polarization would be measured, the corresponding mea- 
surements could be a useful tool for the QGP detection. Such measurements seem 
to be experimentally feasible at RHIC and LHC provided it is supplemented with 
forward detectors. 

Conclusion 

The possibility to reveal a new scattering mode at the LHC is an intriguing one. 
It would significantly change our picture of hadron scattering and lead to bet- 
ter understanding of the non-perturbative region of QCD. Diffraction and related 
processes are very important for studies of collective, coherent phenomena in 
hadronic interactions. Predictions for the experimental observables in these pro- 
cesses presented in this review and their experimental verifications will certainly 
increase the scope of strong interaction studies at the LHC. 
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